Three-Photon Correlations in a Strongly Driven Atom-Cavity System 
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The quantum dynamics of a strongly driven, strongly coupled single-atom-cavity system is studied 
by evaluating time-dependent second- and third-order correlations of the emitted photons. The 
coherent energy exchange, first, between the atom and the cavity mode, and second, between the 
atom-cavity system and the driving laser, is observed. Three-photon detections show an asymmetry 
in time, a consequence of the breakdown of detailed balance. The results are in good agreement with 
theory and are a first step towards the control of a quantum trajectory at larger driving strength. 



Open quantum systems far from thermal equilibrium 
hold great promise for the investigation of fundamental 
physics and the implementation of practical devices 
The versatility of such systems comes from two features: 
the coherent evolution induced by the driving and the dis- 
sipation enabling a transfer of information to an observer. 
These two characteristics affect each other, and the deter- 
ministic evolution is interrupted by unpredictable quan- 
tum jumps 0,0] ■ Monitoring such a quantum trajectory 
is a challenge, in particular when many quantum states 
must be discriminated from each other. A model sys- 
tem in this context is provided by optical cavity quan- 
tum electrodynamics (QED) in the regime of strong light- 
matter coupling, where atomic andphotonic observables 
have been tracked in real time [5|-l8| and controlled by 
means of feedback 0-ll|. However, these experiments 
were performed at low excitation. Stronger driving and, 
hence, faster probing would allow one to track the system 
more closely and explore high-intensity effects like the 
coherent coupling of the system with the drive laser or 
the dynamical polarization of the dressed states [13, [3] ■ 
Moreover, higher excited states containing several pho- 
tons should be discernible by characteristic patterns of 
multiple-photon emissions which can be asymmet- 
ric in time due to the predicted breakdown of detailed 
balance In this Letter we explore such patterns for 
a strongly driven atom-cavity system, when the excita- 
tion rate exceeds the dissipativc rates. 

We consider a system where the atom-cavity coupling 
strength go exceeds the atomic polarization decay rate 7 
and the cavity field decay rate k. The internal dynam- 
ics, described by the Jaynes-Cummings Hamiltonian, is 
complemented by driving with a probe laser of strength 



77 [3l| and dissipation due to spontaneous emission and 
cavity decay. By monitoring the photon stream from the 
cavity, one can evaluate different observables such as the 
average photon number, (a^a), or the average number of 
photon pairs, {a^'^a'^). The first is interesting, e.g., in the 
context of normal-mode spectroscopy [I6, 17 1, while in- 
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FIG. 1: (Color online) (a) A calculated quantum trajectory 
of (a^a) (red) and ) (blue), (b) A zoom reveals coher- 

ent oscillations at different frequencies accompanied by irre- 
versible quantum jumps due to spontaneous emission (V) and 
cavity decay (a). (c) Experimental setup: Single *^Rb atoms 
pass through the cavity mode. Probe light and stabilization 
light are separated using an interference filter (IF). The probe 
light is split by beam-splitters (BS) into four equal parts each 
directed onto a single-photon counting module (SPCM). (d) 
Energy levels of the strongly coupled atom-cavity system with 
the dressed eigenstates |n,±) = {\n, g) ± |n — l,e))/\/2 and 
the coupling (green arrows) induced by the probe laser. 



to undergo complex dynamics as illustrated by the calcu- 
lated quantum trajectory depicted in Fig. [T](a) and (b), 
cf. also [3]. It displays coherent oscillations at different 
frequencies interrupted by sudden quantum jumps due 
to spontaneous emission (V) and cavity decay (a). In 
the following, the time evolution of (a^a) conditioned 
upon such a quantum jump is studied experimentally 
by evaluating the time-dependent second-order intensity 
correlation function. In addition, we introduce measure- 
ments of the third-order intensity correlation function as 



sight into quantum effects can be obtained by regarding 

photon pairs 3, l3- Due to the interplay of the differ- a new tool to study the time-dependence of (a^^^a^) and 



ent dynamical processes, both observables are expected to investigate the dynamics of the system conditioned 
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upon two simultaneous or successive detection events. 
While second-order correlations are naturally symmetric 
in time, third-order correlations enable us to address the 
simple, but experimentally unexplored, question whether 
the emitted photon stream is symmetric in time. 

The experimental setup is depicted in Fig. [l](c) and 



has been described before [llj. A Fabry-Perot resonator 



with a finesse of 195000 and a length of 260 /im, yield- 
ing «;/27r = 1.5 MHz, is tuned to the cycling transition 
from F = 3, mp~3 to F~4, mp = A of the D2-line 
of ®^Rb with 7/27r = 3 MHz. The atom-cavity coupling 
to the TEMqo mode with maximal coupling strength 
f/o/27r = 16 MHz puts the experiment well into the strong 
coupling regime. A laser at 785 nm, detuned by four free- 
spectral ranges, is used to stabilize the cavity length such 
that the bare atom is on resonance with the cavity mode. 
The coupled system is driven by a circularly polarized 
probe laser at 780 nm. In order to excite higher-order 
dressed states, cf. Fig. [l](d), and to increase the signal, 
we use relatively strong probe powers. Since trapping 
atoms is rather difficult in this parameter regime, we per- 
form the measurements with atoms passing through the 
mode, launched via an atomic fountain from underneath 
the cavity. The transit time of an individual atom is 
about 20 /iS. The attractive potential induced by the sta- 
bilization light at 785 nm (AC-Stark shift 5 MHz) guides 
the passing atoms towards regions of strong coupling. 
With a detuning of the probe laser with respect to the 
cavity of Ac/27r = -12 MHz, near-resonant with the lower- 
frequency normal mode, cf. Fig. [T](d), a transient atom 
causes an increase of the probe light transmission. We 
evaluate the recorded photon clicks only in those intervals 
where the transmission is increased by at least a factor of 
1.6 compared to the empty cavity value. In each launch 
about 25 atoms cause such an increase. In the presence 
of an atom, the probability of having a second atom in 
the cavity is less than 3 %. 

First, we consider the time evolution of the average 
photon number shown as a red line in Fig. [T](a) and (b). 
Two characteristic frequencies are visible, a slow oscilla- 
tion with a period around 150 ns and a fast oscillation 
with a period of about 30 ns. To observe these in the ex- 
periment, we evaluate the second-order correlation func- 
tion, (7^^^(r) = {a^ {T)a{T)a) / {a^ a)"^ , of the transmitted 
probe light measuring the conditional evolution of the 
average photon number after the detection of a photon. 

Experimental results are plotted in Fig.[2](a) for differ- 
ent driving strengths (black). Deviations of the asymp- 
totic values from 1 result from the variation of the trans- 
mitted intensity during the passage of an atom through 
the mode and are well understood |2d, 211. Also shown 



are calculations (red), for which details can be found in 
the supplementary information. Since the atomic tran- 
sit happens on a much longer timescale than the internal 
dynamics, we can account for it in the theory by averag- 
ing the correlation function over a proper atomic position 
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FIG. 2: (Color online) (a) 
ent driving strengths (from top 
1.9ft;, 2.7ft;, S.Sk, 3.9ft;, 4.5ft;) in the experiment (black) 
and theory (red), (b) The two oscillation frequencies uj (A) 
and $1 (o) as a function of the driving strength, uj is constant 
with an average tD/27r = 29MHz (dashed blue line) close to 
the maximum vacuum Rabi frequency of 2(7o/27r = 32 MHz. 
Q, increases with rj. Up to a driving strength of 3k (dashed 
red line) it can be described by a two-level model with 
Vl = y/2ri (solid red line). 



distribution, which is the same for all calculations in this 
work. After an additional vertical scaling of the theoreti- 
cal curves by up to 10 % to match the asymptotic values 
of the experiment, we find good agreement between the- 
ory and experiment. 

We identify two different oscillation frequencies in the 
correlation functions. For a quantitative analysis, we fit 
an exponentially damped, oscillating function j32| . The 
obtained frequencies are plotted in Fig. [2](b) as a func- 
tion of the driving strength rj. The faster oscillation fre- 
quency (a) reflects the coherent exchange of energy 
between the cavity mode and the atom, i.e. the vacuum 
Rabi oscillations [22|, |23[ . We find an almost constant fre- 
quency with an average of (D/27r = 29MHz (dashed blue 
line). Due to the atomic motion this is slightly smaller 
than the maximum expected value of 2(7o/27r = 32 MHz. 

The strong driving gives rise to another coherent pro- 
cess, namely the exchange of energy between the atom- 
cavity system and the drive laser. This results in another 
characteristic oscillation frequency (•) which depends 
on the driving strength. As these oscillations are the dy- 
namic manifestation of the supersplitting of the vacuum 
Rabi resonance [2^, we will refer to them as super Rabi 
oscillations. Due to the anharmonic energy-level struc- 
ture of the system, cf. Fig. [1] (d), it behaves at low 
excitation like a driven two level system [3| • Neglecting 
the atom-cavity detuning due to the AC-Stark shift of 
the stabilization laser and the small detuning from the 
normal mode, we expect a Rabi frequency of fi = y/2rj, 
which is plotted as a solid red line in Fig. [2](b). Devia- 
tions from the two-level approximation occur when the 
transition to the second dressed state becomes impor- 
tant. This is expected for a driving strength exceeding 
3k, marked as a vertical dashed line. The reduction of the 
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FIG. 3: (Color online) (a) g'-^\±T,0) (7?=4.5«:). The vacuum 
Rabi oscillations at a frequency of 2go are visible (red arrow) 
as a weak modulation for negative times (red, •) whereas 
the quantum Rabi oscillations at a frequency of 2\/2<?o are 
visible (blue arrow) as a weak modulation for positive times 
(blue, o). (b) A sample quantum trajectory shows the same 
characteristic frequencies for the observables (a^a) (red) and 
(a+=a') (blue). 



oscillation frequency compared to the two-level approx- 
imation at higher powers is in agreement with previous 
studies of a driven anharmonic oscillator (25| . 

Next, we consider the dynamics of another observable, 
namely the probability for the emission of a photon pair, 
{a^'^a'^), blue line in Fig. [T](a) and (b). This is moti- 
vated by the fact that the dynamics of the average photon 
number, (a^a), seems to be dominated by the coherent 
internal dynamics (vacuum Rabi oscillations) and driving 
(super Rabi oscillations) of the first-order dressed states 
only. The quantum Rabi oscillation [27| at a frequency of 
2^/2g of the second-order dressed states are not visible. 
As these states emit photons in pairs IJ] the probability 
to detect a photon pair is more sensitive to the occupa- 
tion of the second-order dressed states than the average 
photon number. This is supported by the corresponding 

(at2a2) 



quantum trajectory. Here, the expectation value 
also undergoes super Rabi oscillations at the same fre- 
quency as (a^a) but with a higher visibility. However, 
the fast oscillations, nicely visible in Fig. [T](b), clearly 
deviate in frequency and visibility from the vacuum Rabi 
oscillations that appear for (at a). 

To confirm this behavior experimentally, we evalu- 
ate the probability to detect a photon pair at a time 
T after a single photon has been observed, correspond- 
ing to the third-order correlation function g^^\T,0) — 
{a^a^{Tfa{Tfa)/{a^a)'^ HI]. For t>0, it measures the 
time dependence of {a^'^a?) conditioned upon the detec- 
tion of a single photon. For r < 0, it measures the time 
dependence of (a^a) conditioned upon the detection of 
a photon pair. A similar time dependence as for the 
second-order correlation function is expected in the lat- 
ter case. 



Experimental and theoretical results are plotted in 
Fig. [3](a) showing good agreement. It is instructive to 
compare these correlation functions directly to a sample 
quantum trajectory shown in Fig.[3](b) where r = is de- 
fined by the detection of a photon. The correlation func- 
tion behaves differently for positive and negative times. 
On very short timescales, a weak modulation at different 
frequencies is clearly visible. For r < 0, a local maxi- 
mum at about 30 ns (red arrow) appears resulting from 
vacuum Rabi oscillations between the normal modes as 
observed in (^'^^(r). The position of the peak matches 
with the oscillation period of the average photon number 
in the quantum trajectory which is also marked by a red 
arrow. For r > 0, we find a shoulder at about 20 ns (blue 
arrow) which was consistently reproduced in other mea- 
surements at different detunings and driving strengths. 
It is in good agreement with the oscillation frequency of 
(flt^a^) in the quantum trajectory. This faster oscillation 
frequency is a consequence of the quantum Rabi oscilla- 
tions of the second-order dressed states with an expected 



period of about 27r/2-\/25o = 22 ns [27|, 

For longer times, the onset of the super Rabi oscil- 
lations is also visible in ^'^'^^(r, 0). While its frequency 
seems to be similar for positive and negative times, the 
amplitude is more pronounced in the first case. Both ob- 
servations are in agreement with our previous statement, 
based on the sample quantum trajectory, that (a^^a^) 
and (at a) oscillate slowly at a similar frequency but with 
different visibility. 

Finally, we investigate how patterns of two successive 
photon detections determine the future trajectory of the 
system. As an example, consider the quantum trajec- 
tory shown in Fig.[T](b), and here the two photons which 
arc emitted at about 650 ns and 680 ns. We ask the 
question whether and how the subsequent time evolu- 
tion of the photon number depends on the separation 
between such two photon detections. To give an answer, 
we evaluate the general third-order correlation function 
9^^\ti,T2) = (atat(Ti)at(ri+r2)a(Ti+T2)a(ri)a)/(ata)3. 
It is proportional to the probability to detect three pho- 
tons, with time separations ti between first and second, 
and T2 between second and third photon. An example 
(77 = 3.9k) is shown in Fig. [4]for experiment (a) and the- 
ory (b). The special case ^('^^(t, 0) discussed previously 
can be found on the vertical (for r>0) and horizontal 
(for T<0) axis of this figure. Our question can now be 
answered by comparing different horizontal cuts through 
the figure. Each cut measures the time evolution of the 
average photon number after the detection of a photon 
pair separated by ti. We see immediately that there is 
a strong dependence on ri, i.e. on the previous mea- 
surement record, if this separation is shorter than the 
coherence time of the system. 

The most interesting feature is the occurrence of a peak 
marked by two dashed lines appearing at ri « 27r/2(7o 
and r2~27r/2r2, i.e. when the separation between the 
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FIG. 4: (Color online) g {ti,T2) in the experiment (a) and 
theory (b) (77 =3.9 k). To reduce noise, it was convoluted 
with a two-dimensional Gaussian filter with a width (FWHM) 
of 7 ns. An intuitive explanation for the pronounced, time- 
asymmetric peak marked by the dashed lines is given in the 
text. 

first two clicks is one period of a vacuum Rabi oscilla- 
tion whereas the separation between the second and the 
third click is half a period of the super Rabi oscillations. 
These timescales suggest that we observe an interplay 
between the vacuum Rabi oscillation and the super Rabi 
oscillation. We give an intuitive explanation in terms of 
the quantum measurement process: The detection of two 
photons separated by risa27r/25o is a signature of the 
normal modes which oscillate at this frequency. Between 
the two detection events the state of the system there- 
fore had a large contribution from the first-order dressed 
states. The detection of the second photon projects these 
states onto the ground state. Subsequently, the super 
Rabi oscillations cause a peak of the excitation after half 
a super Rabi oscillation period, i.e. T2 « 27r/2il. The re- 
verse process at ti sa 27r/2r2 and T2 ~ 27r/2(7o is not partic- 
ularly enhanced: the super Rabi frequency is not charac- 
teristic to any particular set of dressed states. Therefore, 
the detection of two photons separated by half a period 
of the super Rabi frequency docs not change the subse- 
quent time evolution much compared to the detection of 
just an individual photon. This explains the missing of 
a similar peak at ti « 27r/2r2 and T2 ~ 27r/2(7o- 

Both discussed effects, the different dynamics of the 
first- and second-order dressed states as well as the in- 
terplay between vacuum and super Rabi oscillations, give 
rise to pronounced time-asymmetries, i.e. g^^\Ti,T2) 7^ 
g^^H'''2,'''i), demonstrating that the transmitted photon 
stream is asymmetric in time. This effect cannot be ob- 
served for the individual parts of the system, neither the 
empty cavity nor a two-level system in free space. It con- 
firms previous evidence for an asymmetry in studies of 
intensity- field correlation functions Q. The occurrence 
of such time asymmetric fluctuations in the output fields 
can be considered as a direct evidence for the breakdown 
of detailed balance in a driven system far away from ther- 
mal equilibrium [Tsf . 

In conclusion, using second- and third-order intensity 



correlations, we were able to probe the dynamics of the 
normal modes and the second rung of dressed states. All 
relevant dynamical processes - dissipation, driving and 
internal dynamics - have been observed. The next step 
will be to use the information gained by two detection 
events to control the system by means of feedback and 
stabilize its state against fiuctuations. Moreover, higher- 
order correlation functions enable a full characterization 
of the system photon statistics and can therefore be used 
to demonstrate the non-classical nature of the higher- 
order dressed states [29| . 
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Note added: During the preparation of the manuscript 
we became aware of an experiment reporting on super 
Rabi oscillations in a circuit QED experiment (soj . 
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1 Master Equation 



In the rotating wave approximation, the driven Jaynes-Cummings Hamihonian 
reads 

H = hAa(7+(7 + hAcU^a + hg{a) U- + aa^) + hri{a + a)) 

where — uja—^i and Ac ~ uJc — uji are the atom-laser and cavity-laser detun- 
ings, respectively, a (a^) and (t_ {cr+) are the photonic and atomic annihilation 
(creation) operators, g is the coupling strength and rj is the driving strength. 
In the presence of dissipation caused by spontaneous emission and cavity decay, 
the time evolution of the system is governed by a Lindblad master equation [1] : 

p = Cp = — [H, p] + K{2apa) — a) ap — pa) a) + ^{2a-pa+ — cr+a-p — pa+U-) 
in 

where 7 is the atomic polarization decay rate and k is the cavity field decay 
rate. The steady-state density matrix, describing the time-averaged state of the 
system, is defined by 

Cpss = 0. 

After truncating the Hilbert space at some finite photon number (in our case 10) 
the steady state density matrix can be calculated and the time dependent Master 
equation can be solved by diagonalizing it numerically . We use the quantum 
optics toolbox, a Matlab library written by Sze Tan [3], for this purpose. In the 
following, we adapt the formal notation 

p{t) = e^'poiO) 

referring to an explicit solution to the master equation with the initial state pq. 

2 Quantum Regression Theorem 

With this solution to the master equation at hand, one can use the quantum 
regression theorem [31 to evaluate time dependent correlation functions. For the 
second-order correlation function, it states 

{a'^a^ {T)a{T)a) ~ tr {a'ae^'^ [apssO^] } . 

Reading it from right to left, its interpretation is quite straightforward: apssO^ is 
the density matrix of the system after a photon has been annihilated. The term 
e^"^ [apssfl^] designates its time evolution. Except for the normalization, the 
correlation function thus measures the average photon number after the system 
has been perturbed by the detection of a photon. 

The numerical calculation via the quantum optics toolbox simply follows 
the same steps: determine the steady-state density matrix pss, evaluate apssa\ 
propagate it in time by numerically solving the Master equation as discussed 
above, evaluate the expectation value of the average photon number at time r. 
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Similar arguments hold for the third-order correlation functions. Here, the 
quantum regression theorem reads 

{a^ a\Ti)a'^ {ti + T2)a{Ti + T2)a{Ti)a) = tr {a'l'ae^^^ [ae^^^ [apssO.^ } . 

After the detection of the first photon the system is described by the density 
matrix apssQ^- It is propagated in time {e^"^^ [apssO^]) until a second photon 
is detected at ti. The resulting density matrix is propagated again in time 
{e^'^^ [ae^'^^ [apssO^] a^])- The correlation function is then proportional to the 
probability to detect a third photon at time T2 after the detection of the previous 
photon pair. 

3 Position Averaging 

Our previous discussion enables us to calculate the correlation functions for a 
fixed coupling constant g and fixed atom-cavity detuning A^. In practice, both 
parameters depend on the position of the atom and thus vary during the passage 
of the atom through the cavity mode. We write {a'' a'' {T)a{T)a)^ to denote the 
explicit position dependence of the correlation function. On timescales that are 
much shorter than the transit time of an atom through the cavity mode, the 
correlation function that we measure is then given by 

^(2)^^^ ^ / {a''aHT)a{T)a)^p{~^)d^ 
(/ (ata)^p(^)d^) 

where p{~^) is the probability to find an atom at the position 1^ . An analogous 
expression is valid for the third-order correlation function. 

The integration has to be carried out over some finite volume around an 
anti-node of the cavity mode. The size of the volume depends on the trans- 
mission threshold that is used to detect the passing atoms. If the passage of 
the atoms through the cavity mode was unafi^ected by the near resonant probe 
light and the red-detuned stabilization light, one would expect a homogeneous 
distribution independent of position, i.e. pil^) = Po- However, both laser fields 
cause an attractive potential |;5; guiding the atom towards regions of strong cou- 
pling. Therefore, we increased iteratively the probability to find the atom at the 
cavity center by fine tuning the distribution function such that the correlation 
function matches the experimental data. This adjustment was made for one 
parameter set. The same distribution function was then used to reproduce all 
other experimental results shown in the paper using only minor variations of 
the integration volume and a vertical scaling factor as fit parameters. 
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